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ABSTRACT 

We present the first use of the Gemini North laser guide star adaptive optics (LGS AO) system and 
an integral field unit (IFU) to measure the stellar velocity dispersion of the host of a luminous quasar. 
The quasar PGI426+015 [z = 0.086) was observed with the Near-Infrared Integral Field Spectrometer 
(NIFS) on the 8m Gemini North telescope in the H-band as part of the Science Verification phase of 
the new ALTAIR LGS AO system. The NIFS IFU and LGS AO are well suited for host studies of 
luminous quasars because one can achieve a large ratio of host to quasar light. We have measured the 
stellar velocity dispersion of PG1426+015 from 0.1" to I" (0.16 kpc to 1.6 kpc) to be 217 ± 15 kms" 1 
based on high signal-to-noise ratio measurements of Si I, Mg I, and several CO bandheads. This new 
measurement is a factor of four more precise than a previous measurement obtained with long-slit 
spectroscopy and good, natural seeing, yet was obtained with a shorter net integration time. We find 
that PG1426+015 has a velocity dispersion that places it significantly above the Mbh-c* relation of 
quiescent galaxies and lower-luminosity active galactic nuclei with black hole masses estimated from 
reverberation mapping. We discuss several possible explanations for this discrepancy that could be 
addressed with similar observations of a larger sample of luminous quasars. 
Subject headings: galaxies: active — galaxies: kinematics and dynamics 



1. INTRODUCTION 

The link between central supermassive black holes 
and the properties of their host galaxies has been 
firmly established in both qui escent galaxies and ac- 
tive ga lactic nuclei (AGN s; iKormendv fc Richstone 

19951 IFerrarese fc MerrittJ [20001: iGebhardt etal. 

2000allbl: IFerrarese et all 120011: [Graham et all 120011 : 
Marconi fc Huntl 120031 : lHaring fc Rixi I2004D . "Arguably, 
the tightest correlation is the Mbh-c* relation, which 
relates the mass of the central black hole (A/bh) and the 
stellar velocity dispersion of the host spheroid (<j*). In 
active galaxies, the stellar velocity dispersion is typically 
measured using the Ca II triplet stellar absorption 
features at rest wavelengths of 8498, 8542, and 8662A. 
The most direct and broadly applicable method for 
measuring the central black hole mass in active galaxies 
is through reverberation mapping, where the time delay 
between continuum and emission line variations is used 
as a measure of the radius of the broad-line region 
(-Rblr) an d the width of the H/3 emission line (AV) is 
used as a measure of the g as velocity wit hin the BLR 
dBlandford fc McKed H98l iPetersonl Il993h . Then one 
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uses the virial equation to calculate the mass: 



Af, 



BH 



/ 



R BLR (AV) 2 
G 



(1) 



The scale factor / accounts for the unk nown geome- 
try, ki nematics, and inclination of the BLR. fOnken et ahl 
(|2004L hereafter O04) estimated (/) statistically for the 
AGN population under the assumption that the zero- 
point of the A/bh-o* relation is the same for active and 
inactive galaxies. They derive a value of / = 5.5 ± 1.7 
from a sample of 14 active galactic nuclei (AGNs) with 
masses determined by reverberation mapping and stellar 
velocity dispersions determined using the Ca II triplet. 
Using single-epoch spectroscopy, this scale factor is ap- 
plied in black hole mass estimates for large samples of 
AGNs out to high redshift (see McGill et al. 2008 for a 
recent summary). 

The calculation of (/) has unfortunately been lim- 
ited to low-luminosity AGNs for two reasons: dilution 
of stellar features by the AGN continuum and the rel- 
ative scarcity of higher-luminosity AGNs. In particu- 
lar, observations of higher-luminosity objects necessitate 
observations of higher-redshift objects, yet telluric ab- 
sorption mostly limits the utility of the Ca II triplet to 
z < 0.06. Both of these problems can be somewhat cir- 
cumvented with stellar velocity dispersion observations 
using CO bandheads in the near infrared (NIR) H-band. 
This wavelength range corresponds to the maximum ra- 
tio of the stellar e mission from the host galaxy to th e 
quasar continuum (fWright et al.l[l99llElvis et al.lfl99l . 
Additionally, current adaptive optics (AO) systems func- 
tion well in the NIR. This combination facilitates veloc- 
ity dispersion measurements of luminous quasar hosts, 
which are essential to determine if these objects fall on 
the Mbh-c* relation and/or are characterized by the 
same (/) as lower-luminosity AGNs. 
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A suitable sample of higher-luminosity quasars for such 
a study is the r e verbe ration-mapped sample discussed by 
iPeterson et all (|200l . The 16 PG quasars studied in 
this work are on average 40 times more luminous, with 
black holes 10 times more massive, than the AGNs in 
the O04 sample (the O04 and PG quasar samples have 
average luminosities of log(AL A (5100A)/erg s" 1 ) = 43.6 
and 45.2, respectively, and average black hole masses of 
log(M BH /M ) = 7.75 and 8.65, respectively). Work 
has already begun to measure velocity disp ersions in 
this luminosity regime. iDasvra et alj ((2007, hereafter 
D07) investigated the hypothesized evolutionary link be- 
tween PG quasars and ultra-luminous infrared galaxies 
(ULIRGs) with H-band CO bandhead measurements of 
12 PG quasars. These data were obtained using the 
ISAAC long-slit spectrometer (jMoorwood et al.lll998h on 
the 8m Antu unit of the Very Large Telescope under good 
natural seeing conditions. The H-band CO absorption 
features were detected in most of the quasars, including 
four reverberation-mapped quasars. However, the obser- 
vations were still of faint hosts with significant nuclear 
emission contamination. This resulted in rather noisy 
host-galaxy spectra for some objects and velocity dis- 
persions with uncertainties ranging from 18 to 67kms _1 
(12 - 36% error). 

Over the next few years, integral field units (IFUs) 
combined with increasingly reliable AO systems should 
lead to substantial improvements in studies of the hosts 
of luminous quasars. The primary advantage of using 
an IFU is that more of the light in the region of inter- 
est (for example, we probe the central 3" x 3") is dis- 
persed rather than just the light within a normal sin- 
gle slit (e.g., 0.6" x 120" in D07). Consequently, more 
host-galaxy light from near the galaxy's center is gath- 
ered in a single exposure. AO is a further aid because 
the contamination by the intrinsically point-like quasar 
can be confined to the central few pixels of the image, 
thus minimizing the quasar dilution of the stellar absorp- 
tion features. As a demonstration of this approach, we 
present observations of PG1426+015 obtained with the 
Near-Infrared Integral Field Spectrometer (NIFS) dur- 
ing Science Verification Time for the new ALTAIR Laser 
Guide Star (LGS) AO system at Gemini North. We 
specifically chose PG1426+015 because it had a relatively 
uncertain velocity dispersion measurement in the earlier 
study of D07, yet one which placed it suggestively above 
the Mbh-c* relation (along with two of the three other 
luminous reverberation-mapped quasars). 

In the next section we describe our observations, data 
reduction, and velocity dispersion measurement tech- 
nique. We discuss our results and their implications 
in Section [3] and summarize our findings in Section [3J 
Throughout this paper, we adopt Hq = 70kms~ 1 , fl nl = 
0.3, and fitotal = L 

2. OBSERVATIONS, DATA REDUCTION, AND DATA 
ANALYSIS 

The observations of PG1426+015 and the velocity tem- 
plate HIP 75799 (K0III) were carried out at Gemini 
North on 2007 Februa ry 6-7, and April 29 using NIFS 
([McGregor et alJ 12002). Three additional velocity tem- 
plates - HD 84769 (K5III), V* BU CVn (M1III), and 
BD+23 1138 (M5Ia) - were observed on 2008 February 
12, 14, and 15. The image sheer of NIFS divides the 



3" x 3" field of view into 29 spectroscopic slices. NIFS 
was fed by ALTAIR, the AO system at Gemini North, 
using both a natural and laser guide star (|Herriot et al.1 
2000; iBoccas et al.ll2006ft . The quasar was used as the 
natural guide star for tip/tilt corrections. The average 
AO-assisted seeing during our observations was 0.14" in 
the H band. The average uncorrected seeing, as mea- 
sured by ALTAIR at the end of each exposure, was 0.5". 
Typical Strehl ratios for H-band observations using the 
Gemini LGS AO system are about 10%. 

Our observations were obtained with the H grating and 
the JH filter. The spectral resolution is X/5X = 5710 
(FWHM velocity resolution ~ 50kms _1 ) and there are 
1.9 pixels per resolution clement. The H-band wave- 
length coverage from 1.485 to 1.80/im contains many 
strong atomic and molecular features, including Mg I 
1.488/im, Mg I 1.503/mi, CO(3-0) 1.558/mi, CO(4-l) 
1.578/im, Si I 1.589/zm, CO(5-2) 1.598/zm, and CO(6-3) 
1.619/im. 

We obtained a total of 2.08 hours integration on 
PG1426+015 with sequences of 300 second exposures in- 
terspersed with regular sky observations. For telluric cor- 
rection, an A0V star was observed at least once for every 
hour of quasar or velocity template observations, includ- 
ing sky frames. We also obtained the recommended cal- 
ibration frames. 

The data were processed with Gemini IRAF packages 
and the recommended reduction steps given on the NIFS 
website 7 . We only deviated from this prescription to re- 
move the Brackctt absorption features in the A0V stellar 
spectra. For thi s procedure, we use d the xtellcor_general 
task written bv lVacca et al.l (|2003f ). 

Our final host-galaxy spectrum is the difference be- 
tween a 1" (1.6 kpc) and 0.1" radius extraction. We 
empirically chose the outer radius because there was lit- 
tle signal beyond this point and sky subtraction residuals 
became more evident when it was increased. We chose 
the inner radius because it included as much host-galaxy 
light as possible, yet minimized the quasar contribution. 
The LGS system allowed our inner radius to be substan- 
tially smaller than would have been possible under nat- 
ural seeing conditions. 

We measured the stellar velocity dispersion us- 
ing the penalized pix el fitting (pPXF) method of 

iCappellari fc E mscllcm (20Q4|). As is typical in velocity 
dispersion measurements, this method assumes that the 
host-galaxy spectrum is represented by a convolution of 
a stellar template spectrum and the line-of-sight velocity 
distribution (LOSVD). This routine works in pixel space 
and the LOSVD is determined by x 2 minimization. The 
LOSVD is assumed to be in the form of a Gauss-Hermite 
series, where the coefficients are determined simultane- 
ously, but the solution is biased towards a Gaussian. 

3. DISCUSSION 

Two copies of the host-galaxy spectrum of 
PG1426+015 are shown in the top and bottom 
panels of Figure [TJ In the top section, the smooth 
curve shows the spectrum of the M5Ia velocity tem- 
plate, broadened by the best-fit LOSVD. The bottom 
spectrum shows the host spectrum with the broadened 
K5III velocity template. No systematic velocity offsets 

7 http:/ /www. gemini.edu/sciops/instruments/nifs/ 
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Fig. 1. — Two copies of the observed frame host-galaxy spectrum 
of PG1426+0T5, normalized and offset for clarity. The smooth 
curve shows either the M5Ia (top spectrum) or K5III (bottom 
spectrum) velocity template, broadened to fit the host-galaxy ab- 
sorption features using the pPXF method. The gray bands show 
regions of the spectrum that were excluded from the fit for various 
reasons (see text). 

were present between the atomic and molecular features 
in the host spectrum, so we fit all six absorption 
features simultaneously. The shaded areas designate 
spectral regions excluded from our pPXF fit. We 
excluded the 1.658 — 1.682/mi region because it did 
not show absorption features and had moderate telluric 
contamination, the 1.731 — 1.752/im region because 
it was affected by dense sky emission, and the region 
redward of 1.774//m because of increasingly strong 
telluric absorption features. 

The host-galaxy spectrum has an average signal-to- 
noise ratio (SNR) of about 190 per pixel. This is larger 
than the SNR of the long-slit ISAAC spectrum, which 
was obtained with five hours on-source integration. This 
increase in the SNR for a spectrum acquired in less time 
is evidence of the combined advantages of the NIFS IFU 
and the LGS AO system. 

The velocity dispersion resulting from the M5Ia fit is 
224 kms" 1 , with a reduced % 2 of 0.91, while the velocity 
dispersion resulting from the K5III fit is 209 km s -1 , with 
a reduced x 2 of 0.92. The CO(3-0) feature at 1.693/im 
is better fit by the M5Ia template while the Mg I fea- 
ture at 1.632/im is better fit by the K5III template. The 
Mini fit resulted in a velocity dispersion of 207kms _1 , 
with a reduced x 2 of 0.98. Because the M1III fit was 
somewhat worse than the M5Ia and K5III fits, we do not 
include it in Figure [T] and we use the average of only the 
K5III and M5Ia values as our final velocity dispersion. 
The KOIII was not a good match to the host spectrum. 
iGenzel et all (|2001l ) found that K5-M1 giants and super- 
giants are the best templat es for the young stell ar popu- 
lations in ULIRGs (also see Dasyra ct al. 2006a, b) while 
D07 used K0-M0 giants and supergiants as templates for 
the presumably older stellar populations in quasars. The 
fact that later spectral-type templates provide better fits 
to the host spectrum of PG1426+015 may point to a 
young stellar population in the host galaxy. To deter- 
mine the error in individual template fits, we varied the 



velocity dispersion from its best-fit value while holding 
the velocity and higher-order Gauss-Hermite coefficients 
constant. In all three cases (using the K5III, M1III, or 
M5Ia), the location of A% 2 = 1 gave a la error bar of 
~ 10 kms -1 . As a measure of the error due to template 
mismatch, we use the standard deviation of the veloc- 
ity dispersions determined from the three template fits, 
which is also ~ 10 kms -1 . Our final error bar is the 
quadrature sum of these two errors, resulting in a veloc- 
ity dispersion for PG1426+015 of 217 ± 15 kms" 1 . 

This measurement is within an effective physical aper- 
ture of 1.6 kpc. For comparison, the effective radius (r e ) 
for PG1426+015 is 4.86 kpc (based on a Galfit [Peng et 
al. 2002] analysis of Hubble Space Telescope Wide-Field 
Planetary Camera 2 data). Our aperture is intermedi- 
ate between r c and r e /8, two commonly quoted aperture 
sizes in the literature. We consequently choose to use 
our measured value for the remainder of the analysis but 
note that the velocity dispersion would be 208 kms -1 
if we corrected our va lue using the relation derived by 
|j0rgensen et al.l (|1995f ) for E and SO galaxies. This cor- 
rection does not account for the contribution to the ve- 
locity dispersion that we miss by excluding the central 
0.16 kpc; however, this correction is expected to be small. 
For comparison, D07 measured the velocity dispersion 
of PG1426+015 to be 185 ± 67 kms" 1 . This is about 
0.5<7 lower than our measured value, although these mea- 
surements are consistent. Our higher precision is due to 
the higher SNR and broader wavelength coverage of the 
NIFS data. 

Figure [3 shows that PG1426+015 is significantly above 
the Mbh-c* relation for AGNs with masses determined 
by reverberation mapping. The filled squares represent 
Ca II triplet measurements from O04. Th ey have been 
updat ed with velocity dispersion data from lNelson et ail 
(2004) and improved reverberation -bas ed black hole 
masse s from iBentz et al.l (|2006l l2007| ) and iDennev et al.1 
(2006). The open squares represent H-band CO band- 
head measurements of the PG quasars studied in D07. 
The solid line denotes the iTremaine et al.1 |2002) fit to 
the quiescent galaxy Mbh-c * relation. We have not used 
the iFerrarese fe ForJ ([2005) fit to the quiescent galaxy 
Mbh-c* relation because (/) for the AGN population has 
not yet been computed relative to this fit. The velocity 
dispersion presented in this work moves PG1426+015 to 
the position marked by the open star. With our smaller 
error bar, PG1426+015 is now more significantly dis- 
crepant with the Mbh-c* relation. Note that if we com- 
pare the location of PG 1426+015 relative to the steeper 
IFerrarese fe Ford (|2005h Mbh-c* relation fit (and assume 
that (/) is not significantly different from 5.5), it is some- 
what less significantly discrepant. 

D07 found that three of their four reverberation 
mapped PG quasars lie above the Mbh-c* relation and 
list three possible reasons for this. First, the scale 
factor, /, from Equation [T] could be different for PG 
quasars and the lower-luminosity AGNs in the O04 
sample. In general, differences in scale factors could 
arise if different popul a tions have different in clinations 
(e.g.. IWu fc H an 200 11 Uarvis fc McLurd l2006h. accre- 
tion rates (e.g.. ICollin et al. 1 120061 ) . radiation pressure 
(M arcom "et^j]l200l . or some combination of factors. 
Secondly, measurement errors in either cr* or A/bh could 
spuriously drive these objects to lie above the Mbh-c* 
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Fig. 2. — The Mbh~&* relationship for reverberation-mapped 
AGNs. Filled squares represent AGNs with velocity dispersion 
measurements based on the Ca II triplet. These are from O04, 
updated with additional velocity dispersion data and improved 
reverberation-based black hole masses. Open squares show the PG 
quasars from D07, with velocity dispersions measured from the 
CO bandheads. The open star represents our new measurement, 
while the open s quare to its immediate left is the D07 position of 
this object. The[Tremaine et al. (2002) fit to the quiescent galaxy 
Mbh-c* relationship is shown as a solid line. We assume (/) = 5.5. 

relation. D07 note that underestimates of the veloc- 
ity dispersion by 10% are expected because of quasar 
continuum dilution. We found that continuum dilution 
does not significantly affect our measurement by vary- 
ing our inner extraction aperture to include more or less 
quasar continuum. The velocity dispersions resulting 
from these tests were not significantly different from the 
value derived using the optimal inner aperture. In addi- 
tion, preliminary results show that the black hole mass 
in PG2130+099 — one of the D07 quasars that lies above 
the Mbh-c* relation — may be overestimated by a factor 
of several, based on a new measurement of the reverbera- 
tion time lag (Grier et al., in preparation). Thirdly, D07 
show that small number statistics can at least partly ac- 
count for the PG quasars' location above the Mbh-c* 
relation. 

Two final caveats should be considered. First, this 
quasar is an interacting system, with a companion at a 
nuclear separation of 2.7" (4.4 kpc). We would therefore 
expect the velocity dispersion to be larger than predicted 
by the Mbh-c* relation, or opposite to its observed lo- 
cation. Additionally, CO and Ca II triplet velocity dis- 
persions could give inconsistent results (see Dasyra et al. 
2006b for a discussion of this in the context of ULIRGs 
and Silge & Gebhardt 2003 for a discussion in the context 
of early-type quiescent galaxies). If PG1426+015 has an 
atypically young stellar population in its nuclear region, 
as circumstantially suggested by the need for later-type 



spectral templates, these young stars may not trace the 
full velocity dispersion profile of the bulge and cause an 
underestimate of the true stellar velocity dispersion. 

4. SUMMARY 

We obtained H-band observations of PG1426+015 us- 
ing the NIFS IFU on Gemini North and the newly com- 
missioned ALTAIR LGS AO system. The combination 
of using an IFU and LGS AO provides a unique tool for 
studying the hosts of luminous quasars. This is demon- 
strated by the very high SNR spectrum of the quasar 
host galaxy presented in this work as well as our precise 
measurement of the stellar velocity dispersion. We mea- 
sure the stellar velocity dispersion of PG1426+015 to be 
217 ± lSkms- 1 . 

Our measu r ed velocity dispersion and the 
iPeterson etai] ([2004) black hole mass from rever- 
beration mapping places this quasar significantly above 
the Mbh-c* relation for local, quiescent galaxies. We 
have explored a number of reasons why PG1426+015 
might lie above the Mbk-u* relationship. We can 
determine which of these reasons is correct by mea- 
suring velocity dispersions for the remainder of the 
reverberation-mapped PG quasars. It would also be 
valuable to reinvestigate the black hole masses of these 
objects and determine if there are any systematic 
differences between Ca II and CO velocity dispersion 
measurements. When these objectives have been 
completed, we will be able to determine if these higher- 
luminosity quasars lie on the Mbh-c* relation and if the 
scale factor, /, is a function of luminosity. Ultimately 
we hope to include the higher-luminosity sample in the 
calculation of (/) such that future black hole masses 
calculated using this factor will be more reliable. 
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